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Within the resolved Pomeron model of hard diffractive scattering, we compute prompt photon 
production in double-Pomeron-exchange events in proton-proton collisions. Using specific kinemat- 
ical constraints chosen according to the acceptances of the forward proton detectors of experiments 
at the Large Hadron Collider, we provide estimates for inclusive and isolated photon production. 
This is done using the JetPhox program. We find that next-to-leading order corrections to the hard 
process are important and must be included in order to correctly constrain the quark and gluon 
content of the Pomeron from such processes at the LHC. 


I. INTRODUCTION 

A large effort has been devoted to understand the QCD dynamics of hard diffractive events in hadronic collisions, 
since such processes were first observed at HERA mm and at the Tevatron mm more than 20 years ago. Describing 
diffractive processes in QCD had been challenging for decades, but the presence of a large momentum transfer in these 
events brought hope that one could be able to understand them with weak-coupling methods. However, while many 
years of phenomenology point to the existence of a colorless object called the Pomeron, responsible for diffractive 
events when exchanged in the t-channel, understanding hard diffraction in QCD and describing the Pomeron as a 
structure composed of quarks and gluons remains a challenge. 

In the case of deep inelastic scattering (DIS) 7 *p^X, where leptons collide with protons at high energies through 
the exchange of a high-virtuality photon, the situation has reached a satisfactory level. Due to several years of 
experimental efforts at HERA, the diffractive part of the deep inelastic cross-section, which corresponds to about 10% 
of the events, has been measured with good accuracy mm- On the theoretical side, the collinear factorization of the 
DIS cross section also holds for its diffractive component [8] , which allows to separate the short-distance partonic cross 
section computable in perturbation theory from the long-distance dynamics encoded in diffractive parton distribution 
functions (pdfs). 

By contrast, the description of hard diffraction in hadron-hadron collisions still poses great theoretical problems. 
Indeed, Tevatron data provided evidence that even at very large momentum scales, collinear factorization does not 
apply in such cases [9]. In order to estimate hard diffractive cross sections when factorization does not hold, (a modern 
version of) the resolved Pomeron model [10] is being widely used. It makes use of the diffractive pdfs extracted from 
HERA, which give the distribution of quarks and gluons inside the Pomeron depending on the x and kinematical 
variables, while modeling the additional soft interactions that violate factorization. To better test the validity of this 
model, and to better understand the Pomeron structure, it is essential to find sensitive observables to be measured in 
the current colliders experiments. 

One way to constrain quarks and gluons inside the Pomeron is to measure prompt photons in diffractive proton- 
proton (p+p) collisions, as was suggested in [TT]. However, this study relied on leading-order matrix-elements, since 
the Eorward Physics Monte Carlo generator [T2| was used. Subsequent works also relied on LO matrix elements na. 
In this letter, we want to investigate the effects of higher-order corrections, and their impact for a center-of-mass 
energy of 13 TeV at the LHC, and we shall use instead the JetPhox Monte Carlo M to compute the matrix elements 
at leading order (LO) and at next-to-leading order (NLO). 

On the theoretical side, prompt photons refer to high-pt photons created in a hard process, either directly (direct 
photons) or though the fragmentation of a hard parton (fragmentation photons) [15]. On the experimental side, 
inclusive and isolated photons denote prompt photons measured without or with an isolation cut, respectively. These 
are two observables that we shall estimate for double-Pomeron-exchange (DPE) events - meaning diffractive p-fp 
collisions from which both protons escape intact - taking into account the kinematical constraints of the forward 
proton detectors of the CMS-TOTEM Collaborations, or those to be installed by the ATLAS Collaboration in the 
future [T6] . 
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FIG. 1. Leading-order diagrams for prompt photon production in double-Pomeron-exchange events in p+p collisions. Left: 
the annihilation partonic sub-processes are only sensitive to the quark content of the Pomeron. Right: the Compton partonic 
sub-processes are sensitive to the gluon content of the Pomeron as well. 


The plan of the letter is as follows. In section 2, we recall the resolved Pomeron model and its ingredients, explain 
how to obtain diffractive pdfs that take into account the restricted phase space of the outgoing protons, and describe 
how these effective diffractive pdfs are used together with the JetPhox program in order to compute prompt photon 
production in DPE processes in p+p collisions. In section 3, we present our results for a center-of-mass energy of 13 
TeV at the LHC, while analyzing the contributions of Compton, annihilation and fragmentation processes, at LO and 
NLO. Section 4 is devoted to conclusions. 


II. THEORETICAL FORMULATION 
A. Resolved Pomeron model 

The resolved Pomeron model is a long-distance/short-distance collinear factorization framework commonly used to 
calculate hard single-diffract ion and DPE processes. In this work, we focus on DPE prompt photon production at the 
LHC. Related processes (diffractive production of virtual photons and Z bosons) have also been considered previously 
[IIlIIH]. In the case of direct photons, the leading-order diagrams for this process are pictured in Eig. and the cross 
section in the resolved Pomeron model reads: 

= Sdpe E //i/p(^i>^i>/?i>P)/i/p(6,i2,/?2,P) (1) 

where da is the short-distance partonic cross-section, which can be computed order by order in perturbation theory 
(provided the transverse momentum of the photon is sufficiently large), and each factor denotes the diffractive 
parton distribution in a proton. These are non-perturbative objects, however their evolution with the factorization 
scale /i is obtained perturbatively using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi [19] evolution equations. 

In 0 . the variables ^ 1^2 and ti ^2 denote, for each intact proton, their fractional energy loss and the momentum 
squared transferred into the collision, respectively. The convolution is done over the variables I3i^2^ = CiA and 

^2 =^ 2/^2 being the longitudinal momentum fractions of the partons i and j respectively, with respect to the incoming 
protons. However, hard diffractive cross sections in hadronic collisions do not obey such collinear factorization. This 
is due to possible secondary soft interactions between the colliding hadrons which can fill the rapidity gap(s). Eormula 
0 is reminiscent of such a factorization, but it is corrected with the so-called gap survival probability Sdpe which is 
supposed to account for the effects of the soft interactions. Since those happen on much longer time scales compared 
to the hard process, they are modeled by an overall factor, function of the collision energy only. This is part of the 
assumptions that need to be further tested at the LHC. 
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In our computations, we shall use diffractive pdfs extracted from HERA data [20] on diffractive DIS (a process for 
which collinear factorization does hold) by means of an NLO-QCD fit. These are decomposed further into Pomeron 
and Reggeon fluxes /p,R/p and parton distributions /^/p^r: 

= /p/p(C,i)/i/p(/3,M^) + /M/p(C,i)/i/M(/3,M^) . (2) 

The secondary Reggeon contribution is important only at large values of at the edge of the forward proton detector 
acceptance, and therefore we do not take it into account in the following. Measurements at the LHC will allow to 
test the validity of this further factorization of the diffractive proton pdfs into a Pomeron flux and Pomeron pdfs, as 
well as the universality of those Pomeron fluxes and parton distributions. 


B. Effective diffractive pdfs with experimental constraints 

In the following, we assume the intact protons in DPE events to be tagged in the forward proton detectors of the 
CMS-TOTEM Collaborations, or those to be installed by the ATLAS Collaboration in the future [T6|, called AEP 
detectors. The idea is to measure scattered protons at very small angles at the interaction point and to use the LHC 
magnets as a spectrometer to detect and measure them. We use the following acceptances [21]: 

• 0.015 < ^ < 0.15 for ATLAS-AEP 

• 0.0001 < C < 0.17 for TOTEM-CMS . 

Let us now explain how the diffractive pdfs (© are constrained by those detector acceptances. We denote the 
diffractive quark and gluon distributions integrated over t and ^ by and respectively. These 

effective pdfs are obtained from the Pomeron pdfs gp(/d,/i^) and and from the Pomeron flux Let 

us first integrate the latter over the t variable: 

fp{0= dt with fp/p{^,t) = 2ap(t)-i • (3) 

The parameters of Eq.(|^ are the slope of the Pomeron flux Bf> = 5.5^o.7 GeV“^, and Pomeron Regge trajectory 
(ap(t) = (ap(0) + Qfp t with (ap(0) = 1.111 ± 0.007 and o^p = 0.06lo!o6 GeV“^. The boundaries of the t integration are 
^max = i'^p denotes the proton mass) and tmin = ~1 GeV^. The normalization factor Ap is chosen 

such that ^ X dt /p/p(CA) = 1 at ^ = 0.003. 

Next, to obtain the constrained diffractive pdfs, we convolute the Pomeron flux with the Pomeron pdfs while 
imposing a reduction in the phase space of according to the experimental acceptance of the forward detectors: 

|*max(x, .^max) gt 

-r fpiO ■ ( 4 ) 

Jmax(x, .^min) ^ 

Eor the Pomeron pdfs, we make use of the HERA fit B in m- In Eig. 1^ we show the resulting effective diffractive 
pdfs for both the ATLAS-AEP and TOTEM-CMS constraints. These distributions are built in a way to be easily 
incorporated into the LHAPDE library [22] in the grid format. 


C. Computing double-Pomeron-exchange prompt photon production using JetPhox 

JetPhox is a Monte Carlo generator built to compute hadronic cross sections for the process pp yX using the 
collinear factorization framework. Cross sections are calculated as a convolution of short-range matrix elements, com¬ 
puted at LO and NLO, and long-rang^non-perturbative) parton distribution and fragmentation functions. Therefore, 
within the resolved Pomeron model ffl, this program can also be used to compute the cross-sections pp ^ pp^X in 
DPE events. In order to do this, we must substitute the regular pdfs by our effective diffractive pdfs: 

fi/p{x,p) — J didtd(5 S{x- PO = f[/p{x,iP) , (5) 

and multiply the resulting cross sections by the gap survival probability Sdpe- 
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FIG. 2. These plots represent diffractive proton pdfs for three values of with the constraint that the intact proton fall into 
a forward detector. As a result the distributions vanish for > (fmax, and feature a kink when z crosses ^min- Left: diffractive 
gluon distribution with ATLAS-AFP (top) and TOTEM-CMS (bottom) constraints. Right: the sum of the valence 

quarks distribution fq/v with ATLAS-AFP (top) and CMS-TOTEM (bottom) constraints. 


JetPhox produces both inclusive and isolated photons with momentum pt and rapidity y. In case of the inclusive 
cross section, it sums the direct and the fragmentation contribution in the following way: 


da 

dPfdy 


dd^ 

dpjdy 


^ dz 


da°- 

dPtady< 


-{pt/z,y)D2{z,p), 


(6) 


where dd^ is the hard cross section for producing a parton a=(q,q,g) which will then radiate a high-pt photon during 
its fragmentation into a hadron. D2 is the fragmentation function, the z variable is 2 : = p^/pa^ and we have chosen 
the fragmentation scale to be p. In case of isolated photons, an additional criteria is imposed on the hadronic activity 
surrounding the high-p^ photon, as is discussed later. 

In the following, we use JetPhox to compute the direct and the fragmentation contributions in (§, replacing, as 
explained previously, the regular pdfs by the diffractive pdfs extracted above. Technically, this program calls the pdfs 
from the LHAPDF library [22]; we replaced one of those parton distribution sets in grid format by our diffractive 
pdfs constrained with the kinematical cuts. Our choice of factorization scale is p = pt. 


III. NUMERICAL RESULTS 

In this section, we detail the future measurements to be performed at the LHC, in order to test the resolved Pomeron 
model and to constrain the quark and gluon content of the Pomeron, using photon production in DPE processes. We 
use the Monte Carlo program JetPhox (version 1.3.1) to simulate the results, with 2 x 10^ events per channel. 
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FIG. 3. These figures show, for DPE direct photon production, the relative contributions of the Compton and annihilation 
processes at a function of photon pt. Left: at LO, the Compton process represents about 90% of the differential cross section; 
the contribution of the annihilation process is slightly increasing with increasing pt. Right: at NLO, the Compton process 
dominates around 95% of the differential cross section for all the pt range analyzed. 



FIG. 4. These figures show the pt spectrum (left) and the rapidity distribution (right) of DPE inclusive photons (pp ppjX) 
computed by summing the direct and fragmentation contributions for a center of mass energy of 13 TeV, for both ATLAS-AEP 
and TOTEM-CMS detector acceptances. The squares and the triangles represent respectively the LO and NLO calculations, 
the latter giving cross section about 20% greater than the former. 


A. DPE inclusive photons 


In the inclusive mode, there are significant contributions from both direct and fragmentation photons; let us 
first focus on the direct photons. At LO {aemC^s)^ both annihilation processes qq ^ 7 , and Compton processes 
q{q)g q{q)j^ contribute. Going to higher orders opens up additional gluon-initiated sub-processes: gg qqj (NLO) 
and gg ^7 (NNLO) (but the latter contributes only to 1 % of the events). Analyzing the relative contributions 
between the annihilation and Compton channels in DPE events represents a direct way to assess their sensitivity to 
the quark and/or gluon structure of the Pomeron. This is done in Fig. as a function of the photon transverse 
momentum and using ATLAS-AFP acceptance (very similar results are obtained in the TOTEM-CMS case). We 
observe a large dominance of the Compton processes, which could be expected considering the relative magnitude of 
the diffractive gluon and quark distribution shown in Fig. 
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FIG. 5. These figures show, for DPE prompt photon production at NLO, the relative contributions of the direct and fragmen¬ 
tation processes at a function of photon pt. Left: in the inclusive case, the direct and fragmentation contributions are equal 
at Pt —20 GeV, and the relative contribution of direct processes increases with increasing pt. Right: in the isolated case, the 
direct processes dominate; at pt ^ 20 GeV they represent about 75% of the cross section, and that percentage increases with 
increasing pt. 


This means that extracting Pomeron quark distributions from DPE inclusive photon measurements will first require 
that the Pomeron gluon content is well constrained (for instance using DPE dijet measurements [H]). This is even 
more so, since fragmentation photons (which contribute to almost half of the inclusive cross section as we will see 
below) also come mostly from gluon-initiated process. We display in Fig. (left) the differential cross section for 
the production of DPE inclusive photons as a function of the photon pt^ summing all the channels and comparing 
the results at LO and NLO. We show predictions for both ATLAS-AFP and TOTEM-CMS detectors at 13 TeV. 
In Fig. (right), we show the photon rapidity distribution (for pt > 20 GeV), and we note that the difference in 
magnitude between the LO and NLO calculation is about 20%. Obviously, NLO corrections are not negligible, they 
must be taken into account in order to extract correctly the Pomeron structure from future data. Note that due to 
the vanishing of the effective diffractive pdfs for x><frnax, there are no photons produced at very forward or very 
backward rapidities. 

Finally, to compare these cross sections with the future data from the experiments, we note that the gap survival 
probability factor may have to be readjusted. As advocated in several works [231430] . we have assumed Sdpe — 
but the actual value is rather uncertain and must first be measured. The general strategy will be the following. First 
one must perform a global fit of hard diffractive measurements using HERA diffractive DIS data and future LHC 
DPE dijet data, which, if successful, will constraint well the Pomeron gluon content and determine the actual value 
of Sdpe- Then, DPE inclusive and isolated (see below) photon data could also be included in the global fit, which 
would constrain the Pomeron quark/antiquark content, given the increased sensitivity of such processes to those 
distributions. Of course, the level of uncertainty will ultimately depend on how precisely such measurements can be 
performed. 


B. DPE isolated photons 

Using the inclusive photon measurement discussed above in order to constrain the quark content of the Pomeron 
is not optimal, because this observable is contaminated by fragmentation photons, which mostly come from gluon- 
initiated process. In order to suppress the contribution from fragmentation processes, one can use an isolation criteria 
that will disregard the photons that are surrounded by too much hadronic activity. Indeed, generically a direct 
photon will be isolated from a large hadronic activity while a fragmentation photon won’t be, since those are collinear 
emissions from their parent parton. The isolation criteria we use in the following is to require that the hadrons 
measured within a cone of radius R = 0.4 have a maximum transverse energy of 4 GeV. This is one of the options 
available in JetPhox [14], we checked that our conclusions are independent of this particular choice for the isolation 
criteria. 
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FIG. 6. These figures show the pt spectrum (left) and the rapidity distribution (right) of DPE isolated photons for a center 
of mass energy of 13 TeV, for both ATLAS-AFP and TOTEM-CMS detector acceptances. The differential cross sections at 
LO (squares) or NLO (triangles) are obtained by summing the direct and fragmentation contributions while requiring that the 
hadrons measured within a cone of radius 0.4 around the photon have transverse energy no greater than 4 GeV. The NLO 
cross sections are about 50% greater than the LO ones when such isolation criteria is applied. 


Let us now compare the inclusive and isolated photon production. In Fig. we show in both cases the relative con¬ 
tribution of the direct and the fragmentation processes at NLO. The relative contribution of fragmentation processes is 
decreasing with increasing pt^ but it remains always large in the inclusive case: between 20 and 150 GeV, it goes from 
50% to 25%. In the isolated case however, it is clear that the fragmentation contribution is strongly suppressed by 
the isolation criteria, and as the transverse momentum of the photon increases, this contribution eventually becomes 
negligible. 

In Fig.[^ we display our predictions for the pt spectrum (left) and the rapidity distribution (right) of DPE isolated 
photons, for both ATLAS-AFP and TOTEM-CMS detectors at 13 TeV. Comparing the LO and NLO results, we note 
that the NLO cross sections are about 50% greater than the LO ones, which is a much bigger difference than in the 
inclusive case. In order to extract correctly the Pomeron structure from future data, NLO corrections are not even 
more crucial when the isolation criteria is applied. 


IV. CONCLUSION 

Within the resolved Pomeron model 0 . we have analyzed prompt photon production in DPE processes in p+p 
collisions. We have performed our calculations of inclusive and isolated photons using the JetPhox program. This is 
done by substituting the regular pdfs by the diffractive pdfs Q which also take into account the acceptance of the 
forward proton detectors. Then, in order to obtain the DPE cross section, we have multiplied the results by the gap 
survival probability Sdpe^ which we have assumed to be 0.1 for a center of mass energy of 13 TeV. 

Using JetPhox allows us to compute the DPE prompt photon production cross sections with, for the first time, 
next-to-leading order hard matrix elements. Our main result is that the NLO cross sections are larger than the LO 
ones, by about 20% in the inclusive case and 50% in the isolated case. NLO corrections are therefore crucial in such 
processes at the 13 TeV LHC. In addition, we observed that the isolation criteria is necessary in order to suppress 
the contribution of fragmentation photons, radiated by high-pt partons during their fragmentation, and to access in 
a clean way the direct processes of photon production. 

We have also showed that in DPE direct photon production, the Compton partonic sub-processes (Fig [fright) 
clearly dominate over the annihilation ones (Fig left). This is largely explained by the relative magnitude of the 
diffractive gluon and quark distribution shown in Fig. As a consequence, extracting Pomeron quark distributions 
from DPE prompt photon measurements will first require that the Pomeron gluon content is already well constrained, 
which can be done using for instance DPE dijet production m- 







Finally, we have analyzed different possible scenarios to be tested by LHC experiments ATLAS-AFP and TOTEM- 
CMS, and we expect that future data on DPE prompt photon production will provide a quantitative way to test 
the validity of resolved Pomeron model, the factorization of diffractive pdfs into a Pomeron and a Pomeron pdfs, 
as well as to extract the gap survival probability and understand its behavior with increasing energy. We also hope 
that pp^pp^X measurements at 13 TeV will allow to constraint the quark and gluon structure of the Pomeron, as 
pp^^X measurements have helped constrain regular parton distribution functions [3ll[32]. 
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